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Willis and Kerswell Reply: In our letter [1] it was 
reported that in pipe flow the median time r for relami- 
narisation of localised turbulent disturbances closely fol- 
lows the scaling r ~ l/(Re c — Re). This conclusion was 
based on data from collections of 40 to 60 independent 
simulations at each of six different Reynolds numbers, 
Re. In the previous comment, Hof et al. estimate r dif- 
ferently for the point at lowest Re. Although this point 
is the most uncertain, it forms the basis for their asser- 
tion that the data might then fit an exponential scaling 
r ~ exp(A Re), supporting [2], for some constant A. The 
most certain point (at largest Re) does not fit their con- 
clusion and is rejected. We clarify why their argument for 
rejecting this point is flawed. The median r is estimated 
from the distribution of observations, and it is shown that 
the correct part of the distribution is used. The data is 
sufficiently well determined to show that the exponential 
scaling cannot be fit to the data over this range of Re, 
whereas the r ~ l/(Re c — Re) fit is excellent, indicating 
critical behaviour and supporting experiments [3]. 

The data median was used in our original analysis as 
the estimator for the parameter of the exponential distri- 
bution. Estimating via the log-probability plot as in [2[, 
the range of possible lines that can be drawn does not at 
all reflect the true error in the mean. In the following we 
estimate confidence intervals by using the bootstrapping 
method Q applied to the maximum likelihood estimator 
of 1/t. Any truncated data is integrated in a logical fash- 
ion into this resampling method by further resampling. 
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FIG. 1: Data with 95% confidence intervals. Inset: median 
and 95% bracketing lines vs. number of observations removed. 

The data are affected by a transient time to for adjust- 
ment from the initial condition at Re = 1900 to the oper- 
ating Re. This is of around 50 D/U in the worst case. In 
the inset to Fig. [T] we calculate 1/r with 95% confidence 
intervals as a function of the number of smallest observa- 
tions ignored. In the original analysis a maximum of 10% 
of the data was removed and for Re = 1580 gives the esti- 


mate indicated by the first solid bar. Cutting more leads 
to a rise in 1/r but the uncertainty becomes very large. 
A possible time as suggested by Hof et al. is r = 4.9 D/U 
(second higher bar, and red x on Re plot), but it is an 
order of magnitude smaller than to. As r <C to so few ob- 
servations survive past the transient that nearly 20 of 40 
observations must be cut. Aside from statistical uncer- 
tainty, this short time is far from the asymptotic r — > oo 
and it is questionable that a self-sustaining mechanism is 
really involved. The exponential fit of Hof et al, however, 
relies heavily on this data point. 

The largest step in Re is from 1580 to 1700. Here the 
adjustment in Re from the initial condition is smaller, the 
transient time is shorter, the typical observation time is 
longer and the 'tail' exponential distribution is clearly 
seen. It is this part of the distribution which has been 
fit - after the first few points the estimate for 1/r (solid 
bar) is within the confidence interval even when large 
proportions of the data are ignored. 

At Re = 1860 Hof et al. state that the data is in- 
conclusive because the probability P(T) is only shown 
for T < 1000 D/U [5]. A maximum observation time is 
not a limitation provided that the probability distribu- 
tion observed with in this time reflects that beyond it. 
The data fit a straight line on the log-plot and therefore 
very well the exponential distribution, i.e. the process is 
memory less. For the behaviour to change to a second 
memoryless process beyond T — 1000 D/U is contradic- 
tory, requiring some property of the flow to be preserved 
to such long times. 

In Fig. [T] the best-fit exponential and inverse scalings 
are fit to data at the 5 largest Re. It is evident that most 
points are significantly different at the 5% level from the 
exponential scaling over this range of Re, no matter how 
it is fit. In contrast l/(Re c — Re) with simple power — 1 
is a very natural fit. 

We thank David Leslie for advice on the bootstrapping 
method. 
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